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Abstract 

Droplet solitons are a strongly nonlinear, inherently dynamic struc¬ 
ture in the magnetization of ferromagnets, balancing dispersion (ex¬ 
change energy) with focusing nonlinearity (strong perpendicular anisotropy). 
Large droplet solitons have the approximate form of a circular do¬ 
main wall sustained by precession and, in contrast to single mag¬ 
netic vortices, are predicted to propagate in an extended, homoge¬ 
neous magnetic medium. In this work, multiscale perturbation theory 
is utilized to develop an analytical framework for investigating the 
impact of additional physical effects on the behavior of a propagat¬ 
ing droplet. After first developing soliton perturbation theory in the 
general context of Hamiltonian systems, a number of physical phe¬ 
nomena of current interest are investigated. These include droplet- 
droplet and droplet-boundary interactions, spatial magnetic field in¬ 
homogeneities, spin transfer torque induced forcing in a nanocontact 
device, and damping. Their combined effects demonstrate the funda¬ 
mental mechanisms for a previously observed droplet drift instability 
and under what conditions a slowly propagating droplet can be sup¬ 
ported by the nanocontact, important considerations for applications. 

This framework emphasizes the particle-like dynamics of the droplet, 
providing analytically tractable and practical predictions for modern 
experimental configurations. 


1 Introduction 

The ability to excite, probe, and control magnetic media at the nanometer 
scale has enabled new applications in spintronics and magnonics as well as 
the exploration of new physics. Of particular, recent interest is the genera¬ 
tion of coherent structures within the magnetization, a vector held describ¬ 
ing the magnetic dipole moments of a magnetic material. Coherent struc¬ 
tures observed experimentally include those with nontrivial topology such 


1 



as magnetic vortices Uhlir et al. 2013|, 

Pulecio et al. |2014 and skyrmions 

Nagaosa and Tokura 

[2013 , 

Fert et al. 

2013 , Yu et al. 2014 as well as 

nontopological droplet solitons Mohseni et al. 

[2013 2014 , |Chung et al. 


|2014 , Macia et al. [2014 . The magnetic droplet soliton (droplet hereafter) 
is a coherently processing, nanometer-scale localized wave structure exhibit¬ 
ing strongly nonlinear effects Kosevich et al. |l990 


Its first observation 


Mohseni et al. |2013| was enabled by a nano-contact spin-torque oscillator 


(NC-STO) device, which provides the necessary forcing to oppose magnetic 


damping, hence has been termed a dissipative droplet soliton Hoefer et al. 
[2010 


The droplet is theoretically understood as a solution of a conservative 
Landau-Lifshitz partial differential equation modeling spatio-temporal dy¬ 
namics in an ultra thin, two-dimensional magnetic him with uniaxial, per¬ 
pendicular anisotropy. When the stationary droplet’s processional frequency 
is close to the local, magnetic held induced (Zeeman) frequency, it resembles 
a circular domain wall of large radius, almost reversed at its core. These 
large droplets can propagate, coinciding with a superimposed wave struc¬ 
ture to the otherwise spatially homogeneous phase. An example is shown in 
Fig. [T] 



Figure 1: (a) Contour (color is m z ) and vector held, (m x ,m y ), plot of an 
approximate, propagating magnetic droplet, described in Sec. [3| moving to 
the right with speed V = 0.02 and processing with frequency u> = 0.1. (b) 
The magnetization evaluated on the line y = 0 and superimposed upon the 
m z prohle. There is no phase winding across this nontopological structure. 


There are a rich variety of physical mechanisms that can change the ori¬ 
entation of the local dipole moment in a ferromagnet and hence inhuence the 
particle-like behavior of a droplet. An example already mentioned is that 
of the NC-STO and damping. This motivated us to develop a droplet soli- 
ton perturbation theory that allows for the analysis of droplet dynamics in 


the presence of a large number of physical perturbations Bookman and Hoe¬ 


fer] [20131. This theory describes soliton dynamics via a finite dimensional 


dynamical system representing the adiabatic evolution of the droplet’s four 
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parameters (center, precessional frequency, and phase) resulting from per¬ 
turbations. 

While the theory provided fundamental explanations of droplet physics, 
it was limited to almost stationary droplets, i.e., droplets of negligible mo¬ 
mentum, which manifested as a constraint equation on the dynamics. Fur¬ 
thermore, important physical effects such as droplet acceleration due to a 

p98 


magnetic field gradient Kosevich et al. 


Babich and Kosevich 


Hoefer et al. |2012| and droplet interactions Piette and Zakrzewski 
Maiden et al. |2014j were excluded from the theory. 


2001 


1998 


Soliton perturbation theory has been successfully used to describe dy¬ 
namics in many physical systems Kivshar and Malomed |1989| , Sanchez and 


R. 1998 , notably solitons in Nonlinear Schrodinger (NLS) type equations 


modeling, for example, optical fibers Yang 2010 and Bose-Einstein conden¬ 
sates Kevrekidis et al. 2008 . The central idea is to project the perturbed 


PDE dynamics onto the unperturbed soliton solution manifold, allowing for 
adiabatic evolution of the soliton’s parameters. The resulting modulation 


equations can be obtained in different ways, of which multiple scales per- 

turbation theory 

Keener and McLaughlin 1977 or perturbed conservation 

laws Kivshar and Malomed 

[1989 are perhaps the most common. Both ap- 


proaches have been shown to be equivalent in some specific applications, 
see, e.g., Ablowitz et al. [2009], but the conservation law approach is power¬ 
ful in its simplicity. However, it can be unclear which balance laws to use, 
especially when higher order information is sought. The governing Landau- 
Lifshitz equation for magnetization dynamics is a strongly nonlinear, vecto¬ 
rial generalization of the NLS equation that lacks Galilean invariance. There 
exist two-dimensional moving droplet solutions Piette and Zakrzewski |1998| , 
Hoefer and Sommacal [2012 described by six independent parameters. As 


such, droplet perturbation theory is significantly more complex than its NLS 
soliton counterpart, so a structured approach is desirable. 

The earliest, most general formulation of soliton perturbation theory we 
have found in the literature is given in Keener and McLaughlin |1977 , where 
multiple scales perturbation theory and the Green’s function for the lin¬ 
earized operator are utilized. An alternative, rigorous approach was devel¬ 


oped for the modulational stability of solitons in NLS equations in Wein¬ 


stein |l985] based on multiple scales and the generalized nullspace of the 
linearized operator. Motivated by the challenges associated with magnetic 
droplet soliton perturbation theory, we revisit the general approach for per¬ 
turbed Hamiltonian systems, utilizing multiple scales and the generalized 
nullspace formulation. In Section [2j we provide necessary conditions and 
expressions for the resulting modulation equations. Our approach preserves 
much of the generality of the Green’s function approach with the comparative 
accessibility of Weinstein’s approach to NLS. Additionally, we demonstrate 
that these methods do generalize to soliton perturbation in higher dimen¬ 
sion, beyond (1+1)D, and how they can be used to determine the evolution 


3 







































































of large numbers of parameters for a single soliton. 

Two-dimensional moving droplet solutions have been computed |Piette 
and Zakrzewski 1998 , Hoefer and Sommacal |2012 and studied asymptoti¬ 
cally in the weakly nonlinear regime Ivanov et al. 20011. They can be accel¬ 
erated and controlled by an inhomogeneous, external magnetic field Hoefer 
|2012 | and exhibit novel interaction properties Maiden et al. 2014 


et al. 


Droplets can also experience a drift instability in NC-STOs, whose origin is 
not well-understood Hoefer et al. 2010 . Because droplets are relatively new 
physical features of nanomagnetic systems, they hold potential for applica¬ 
tions such as spintronic information storage and transfer or probing material 
properties. Moreover, their fundamental physics are not well understood. A 
more general theory to describe the motion of droplets in realistic physical 
systems is desirable. 

Toward this end, we here present an analytical framework for investi¬ 
gating the impact of a large class of physical effects on a six parameter 
propagating droplet soliton. In Section [3j we derive an approximate solution 
for the propagating droplet in the close to Zeeman frequency (large mass) 
and small velocity (order one momentum) regime, which greatly reduces the 
complexity of the asymptotic theory developed. This allows for explicit ana¬ 
lytical results, a powerful feature of soliton perturbation theory considering 
the droplet’s strongly nonlinear qualities. The main result of this paper, 
the modulation equations, is presented in Section [IJ utilizing the general 
soliton perturbation theory formulation from Section [2j The versatility of 
this framework is subsequently discussed in Sections [5] and [6] through an 
investigation of a series of perturbations of current physical interest. In par¬ 
ticular, we analytically demonstrate a mechanism leading to the NC-STO 
droplet drift instability and why a droplet can be attracted or repelled by a 
ferromagnetic boundary. 

The model under study here is the perturbed Landau-Lifshitz equation 
for the magnetization, m = [m x , m y , m x \, of a two-dimensional ferromag¬ 
netic film in nondimensional form 


—— = —m x h e g + ep, m : M 2 x (0, oo) - 
h e g = V 2 m + (ho + m~) z, lim m = z, 

|x |—>OQ 


(1) 


where p is a perturbation that preserves the magnetization’s length (p • m = 
0) and 0 < e <C 1 is a small parameter encoding the strength of the perturba¬ 
tion. The perturbation can depend explicitly on time so long as the variation 
is slow. That is, the perturbation depends only on the slow temporal coor¬ 
dinate T = et. The perpendicular, external magnetic field hoi = ho(X.,T)z 
can be large and slowly varying in space (X = ex) and time. Further gener¬ 
alizations to rapidly varying perturbations are possible. A necessary require¬ 
ment for the existence of droplet solitons is strong perpendicular anisotropy, 
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encoded in the orientation of the m z z term in the effective field h e fj. See 
Bookman and Hoefer 12013 for the derivation and nondimensionalization of 
this model. 


2 Adiabatic Dynamics for Hamiltonian Systems 


The main aim of this paper is to determine the evolution of parameters of 
the droplet soliton in response to a fairly general class of perturbations. The 
equations determining the time-evolution of parameters will be referred to 
as modulation equations. Before performing that analysis directly, we first 
consider the more general context of perturbed Hamiltonian systems. Note 
that the Landau-Lifshitz equation is a Hamiltonian system, with canonically 
conjugate variables cos(0), <f>. That is, the Landau-Lifshitz system may be 
written as and = — g c '^ q , where the right hand sides are 

expressed in terms of variational derivatives of the energy £, defined later 
in eq. & In such systems, it is more convenient to execute perturbation 
theory in the Hamiltonian variables. It is also possible to write down mod¬ 
ulation equations for solitons with parameters related in a particular way. 
This statement will be made more precise later in this section. 

The basic procedure is to allow the parameters to vary on a time scale 
proportional to the strength of the perturbation, e. By allowing the param¬ 
eters to vary in this way, additional degrees of freedom are introduced which 
can be used to resolve the difficulties arising from singular perturbations. 
Expanding about the soliton solution in an asymptotic series, one obtains a 
linear problem at order e. In general, this linear equation will not admit solu¬ 
tions bounded in time. However, as utilized in Weinstein 11985 , a solvability 
condition exists whereby bounded solutions are assured, guaranteeing that 
the linear problem at order e does not break the asymptotic ordering. Im¬ 
posing these conditions leads to the modulation equations. This procedure is 
equivalent to projecting the solution of the perturbed model onto the family 
of solitons. While one might wish to then solve the linear equation at order 
e to obtain a further correction, we will not do that in this work. As will be 
demonstrated by the examples in Sec. [5j quite satisfactory predictions can 
be made considering only the leading order dynamics. 

A Hamiltonian system requires a real inner product space, X; a nonlinear 
functional, H : X —>• M; and a skew adjoint operator J : X —> X. We will 
use the notation •} for the inner product on X. The standard form for a 
Hamiltonian system is 

^ = JVH(z) (2) 

where 2 £ I is referred to as the state variable. We refer to H as the 
Hamiltonian, which is often assigned the physical meaning of energy since 
it is automatically a conserved quantity of such a system. In this context, 
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by Vi? we mean the first variation of this nonlinear functional and by AH 
we mean the second variation (both taken with respect to the state variable, 
z). We consider here Hamiltonians which depend explicitly upon additional 
parameters, q E M m (m is the number of such parameters). Such parameters 
may arise due to a change of coordinates, such as to a comoving reference 
frame. For the examples which arise in this work, the parameters q arise 
from just such a transformation, so we will typically refer to these parameters 
as “frequencies". 

We assume here that (|2]) admits a solitary wave solution, u. 

0 = JVF(«,q). (3) 


If H depends on q, naturally u will depend on q as well. Typically, the 
parameters q do not provide a full parameterization of the solitary wave 
manifold due to underlying symmetries in the equation such as translation 
invariance. Accordingly, we will allow u to depend on a separate set of 
parameters r E M s (s is the number of such parameters). For reasons that 
will become clear in later examples we refer to these parameters as “phases”. 
Therefore, the solitary wave can be written u = u(x; q, r). Often times, the 
Hamiltonian system ([2]) admitting solitary wave solutions ([3j is idealized, 
neglecting important physical effects. While some perturbations may give 
rise to a different Hamiltonian system, in general such effects do not need 
to preserve the Hamiltonian structure. We will treat both cases the same 
by introducing a small perturbation into the equation itself. The perturbed 
model is 

= JVi?(z,q)+eP (4) 

where 0 < e <C 1 and P is a perturbation. The parameters q, r are allowed to 
vary on a slow time scale, T = et. We restrict to perturbations which depend 
explicitly on time only through this slow time variable, T. In this case, 
ordinary differential equations governing the evolution of these parameters 
can be determined according to the following theorem. 


Theorem 1. Given the perturbed Hamiltonian system @>. If 

1. The solitary wave solution, u, exists for the unperturbed system Q, 
e = 0, and is independent oft. 

2. J is invertible. 


3. AH 


is self-adjoint for all admissible q. 


4- V 1 < k < m, 3 1 < j < s such that VH(z, q) 


_ _ t—1 du 

~ J dr, 
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then letting v = [r, q] T E M s+m , the modulation equations are 


( s-\-m 

E 

V 2=1 




5m 7 5m 7 - / dT 


J-'P. 


5m 
’ dvj 


( 5 ) 


Equation [5] is consistent with previous general results when applied to 
Hamiltonian systems Keener and McLaughlin |1977 . The assumptions of 


Theorem [l] may seem restrictive at first, but these conditions are frequently 
met in physical systems of interest. In all systems under consideration here 
there does exist a solitary wave solution. These solutions generically depend 
on time, but for the case of a single solitary wave solution, transforming 
to the reference frame moving, rotating, and/or processing with the solitary 
wave can eliminate this explicit dependence on time. Such a transformation 
will introduce parameters in q and alter the Hamiltonian but leaves the 
Hamiltonian structure intact. 

The second does offer a restriction. For instance, in the Korteweg-de- 
Vries equation, J does not admit a bounded inverse and correspondingly the 
modulation equations require additional considerations Ablowitz and Segur] 
[1981 . Nevertheless, formal calculations are possible and J is frequently 


invertible for Hamiltonian systems (as it is, e.g., for NLS and the Landua- 
Lifshitz equation). 

With appropriate restrictions on the Hamiltonian, the third assumption 
always holds. The self-adjoint property of the second variation essentially 
follows from the same calculation which proves the equality of mixed partial 
derivatives in finite-dimensional calculus. More care needs to be taken in the 
corresponding calculation on function spaces, but the Hamiltonians derived 
in physically relevant systems typically are well enough behaved. 

The fourth assumption is restrictive and may seem obscure. However, the 
parameters of the soliton are often speeds or frequencies. These parameters 
are typically linked to initial positions or initial phase values so that q and r 
have the same length (s = m). In such cases, the dependence of the soliton 
on the parameters in the laboratory frame will be in the form r + tq. From 
this temporal dependence, the relations in Assumption (iv) follow directly. 


2.1 Derivation of Equation ([5]) 

Theorem [l] relies on the following lemma. 

Lemma 1. Let X be a Hilbert space. Let A be a linear operator mapping X to 
itself. Let f E X. Let A t be the adjoint of A, i.e. the unique linear operator 
satisfying (A^x,y) = (x,Ay) for all x, y E X. Define T : [0, oo) —> X as the 
solution of the initial value problem 


( § = AT + f 
T(0) = T 0 E X. 


( 6 ) 
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Let [ 1 -1 = 0 and = /Xj_i /or 0 < z < N, where N denotes the highest 
integer such that (A^) N has nontrivial kernel. Then T(t) will not be bounded 
in time unless (/ij_i, Tq) + (//j, /) =0 /or 0 < i < N. 


m 


This lemma is a minor generalization of the solvability condition proven 
Weinstein 11985 . There are a few key limitations which may not be clear 
upon first reading the statement of the lemma itself. First, A and / are 
assumed to be independent of time t. Second, all assumptions of smoothness 
of / are bound up in the choice of X which is problem specific. In the context 
of Hamiltonian systems, X is given and the required smoothness of / is clear. 
In our intended application, Eq. ([6j arises from a linearization of a nonlinear 
problem about a given state. In this case, A and / are given, but not X. In 
order that Lemma [l] applies, there must exist an X which makes A and / 
compatible, and it will be in that sense which T(t) remains bounded in time. 
From here on out, we assume sufficient smoothness in our perturbation such 
that a Hilbert space is naturally chosen. For the perturbations we investigate 
in Section [5j this is the case. Finally, the details of defining the adjoint of 
the unbounded operator A are not considered here but can be handled in a 
standard manner; see, e.g., Weinstein |1985| . 

The proof of Theorem [l] proceeds by substituting the ansatz 


^ = u(r(T), q(T)) + eui(x, t,T) + O (e 2 ) 


(7) 


into Q. Expanding in powers of e, the first order equation becomes 

< 8 > 


Note that Eq. Q is of the form in Lemma [I] (A = JAH(u,q), f = P — 
friST ~ or der that the expansion in ([7]) remain asymptotically 

ordered, it is necessary that u±(x,t,T) remain 0{ 1) for sufficiently long 
times. Lemma [I] thus gives a condition that must be satisfied. It remains 
to characterize the generalized nullspace of (JAH(u, q))7 Note that since 
AH(u,q) is self-adjoint, (JAH(u, q))t = — AH{u, q)J. 

Differentiating ([3| with respect to the parameter Tj for 1 < j < s and 
applying J~ 1 to the result yields AH(u,q)^: = 0. It follows that J 

is in the kernel of (JAH(u, q))t for all j. Differentiating ([3]) with respect to 
the parameter qk for 1 < k < m yields 




= 0 . 


(9) 


Utilizing assumption (iv), we can replace the second term in (J9| so that there 
is some j with 


AH(u, q) J J 


T -1 9U 

dqk 


= J 


_ 1 du 
drn 


( 10 ) 
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Hence, £ ker(A H(u. q) J) 2 and therefore in the generalized nullspace. 

These two sets of vectors do not necessarily characterize the full generalized 
nullspace; however, these offer a sufficient number of constraints to uniquely 
determine the modulation system. Requiring that f = P — — fq 5 r^e 

orthogonal to J _1 J^~ and J~ 1 §^ L yields equations (|5]). The modes 
and r: may not give rise to a complete characterization of the nullspace. 
As a result, Eqs. ([5]) are only a necessary but not sufficient condition to 
prevent secular growth. 


3 Approximate Propagating Droplet 


Now that we have the general formulation of the soliton modulation equa¬ 
tions for perturbed Hamiltonian systems in eq. (|5|, we would like to apply 
them to the magnetic droplet soliton solution of eq. ([!]). For this, we will 
need to compute derivatives of the droplet with respect to its parameters 
as well as associated inner products. This could be performed numerically 
with a “database” of droplet solutions as in Hoefer et al. 12012 . Here, we 


obtain an explicit, analytical formulation of the modulation equations in the 
strongly nonlinear, moving droplet regime. But before we can determine the 
modulation equations, we need an explicit representation of the propagating 
droplet itself. In this section, we derive an approximate solution to eq. ([!]) 
when e = 0, a restriction we maintain for the remainder of this section. The 
solution describes a slowly moving droplet with frequency just above the 
Zeeman frequency. A droplet solution can be characterized by six param¬ 
eters: its precession frequency uj above the Zeeman frequency ho in these 
non-dimensional units, propagation velocity V = [V x ,V y ] T , initial phase $0) 
and the coordinates of the droplet center ^ = [£ x , £, y ] T = Vi + xo- 

Approximate droplet solutions have been found in two regimes: (i) 0 < 
1—a; — |V| 2 /4 <C 1, near the linear (spin-wave) band edge correspond¬ 
ing to propagating, weakly nonlinear droplets approximated by the NLS 

20011 elocity 


Townes soliton Ivanov et al. 


corresponding to stationary, strongly nonlinear droplets approximated by 
a circular domain wall Kosevich et al. |1986 , Ivanov and Stephanovichj 
|l98 9|. We will focus here on large amplitude propagating solitons where 


the magnetization is nearly reversed because experiments operate in this 
regime. Note, however, that the weakly nonlinear regime could also be stud¬ 
ied. The defining equation for the droplet can be formulated as a bound¬ 
ary value problem by expressing the magnetization in spherical variables 
m = [sin(@) cos(<I>), sin(0) sin(<I>), cos(©)] in the frame moving and precess- 
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ing with the soliton 0 —> 0(x — £), <h —> 4>o + (ho + u)t + <h(x — £): 

— sin(0)V • V0 
, sin(0)(w- V- V$) 
lim V4? 

„ |x|-K» 

This problem can be further simplified by exploiting the invariance of 
Eq. © under rotation of the domain to align the x-axis with the propaga¬ 
tion direction. In this coordinate system, V = Ex. Adding the assumptions 
of small frequency and propagation speed, a simple correction to the known, 
approximate stationary droplet can be found (see Appendix 0- 


V • (sin 2 0V$)) 

—V 2 0 + \ sin(20)(l + |V$| 2 ) 
— y, lim 0 = 0. 


( 11 ) 


0 = cos 1 



<h = d>0 + (ho + u)t 


V 




cos(<£>) + O 



( 12 ) 

(13) 


Above, (p, tp) are polar variables for the plane, whose origin is centered on the 
droplet. That is p = — £ x ) 2 + (y — £ y ) 2 and tp = arctan ^^ . See 

Fig-0 for a visualization of this approximate solution. This approximation 
is valid so long as 

0<|E|<cj, 0 < w < 1. (14) 


As for the stationary case, the propagating droplet can be viewed as a pre- 
cessing, circular domain wall with a radius that is the inverse of the fre¬ 
quency. The new term — V cos(</?)/co 2 reveals the deviation of the propagat¬ 
ing droplet’s phase from spatial uniformity. While the relations in (14) may, 
at first, seem overly restrictive, we will show that important and practical 
information about propagating droplets can be obtained in this regime. This 
approximate solution offers both an error estimate and is amenable to fur¬ 
ther analysis in the context of the perturbed Landau-Lifshitz equation Q. 
Furthermore, it provides a significant improvement over the approximate 
droplets used in past numerical experiments Piette and Zakrzewski |1998| 
when the asymptotic relations (14) hold. 

Another important property of Eq. ([!]) in the e = 0 case is that it admits 
conserved quantities, including the total spin 


Af 


m z )dx, 


(15) 


the momentum 

V 


rriyS7m x - m x X7m y 

1 + m z X ’ 


(16) 


10 














and the total energy, 


£ = 


1 


(jVm| 2 + (1 — m 2 )^ dx + 


(17) 


These quantities are not independent for the droplet itself, since the droplet is 
the energy minimizing solution constrained by the total spin and momentum. 
Utilizing the approximate form (12), (13) for the droplet, a map can be 


constructed between its parameters and the conserved quantities. Evaluating 
the integrals in Eqs. at the approximate droplet, we obtain 


, r _ 27r 
A w 2 ’ 

V = 

ui 6 

£ = u? (l v l 2 + 4u;2 + 


(18) 

(19) 

( 20 ) 


where higher order terms in ui and |V| have been neglected. These formulae 
extend the predictions for stationary droplets, see, e.g., Kosevich et al. 19901, 
and offer an analogy to classical particle dynamics. Rewriting £ in terms of 
the other conserved quantities, we obtain 


s = Vto( l lv 


2 A7 


3 + M* 


1 

+ • 


( 21 ) 


By analogy to classical systems, we can interpret \/27r |'7-’| 2 /2AT'z as the 
kinetic energy of the droplet, \f2itN 2 as the droplet’s potential energy due 
to precession, and h^N / 2 as the Zeeman energy of the droplet with the net 
dipole moment J\T. Inspection of the kinetic energy term shows that 


m eS = 


M z/2 _ 2vr 

V2n 


LU 


( 22 ) 


serves as the effective mass for the droplet. Therefore, the 0<w<l regime 
corresponds to droplets with large mass. This is a natural interpretation 
since it is the precession of the droplet which determines its size and prevents 


the structure from collapsing in on itself. On the other hand, eq. (19) implies 
that the slowly propagating |V| <C uj regime supports droplets with up to 
\V\ = O (^) momenta. We will return to this observation of an effective mass 


for the droplet in Section ^3.1 where we consider the dynamical equations 
induced by spatial inhomogeneity in the external magnetic field. 

One description of the magnetic droplet is as a bound state of magnons 
(19901. It is then natural to interpret the potential energy 


Kosevich et al. 


as the energy released by decay into these constituent “subatomic 


particles”. The expressions (18) and (19) can also be utilized to verify the 
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Vakhitov-Kolokolov soliton stability criteria Vakhitov and Kolokolov Jl973|, 


Grillakis et al. |1990 for a propagating droplet (see Hoefer and Sommacal 
2012] ), namely that A f u = —47 t/lo 3 < 0 and A/^Vy • 'P — VyAA • P u = 


— 8iT~/uj e < 0, as required. 

For the remainder of this work, we will use the approximate droplet in 


eqs. (12), (13). 


4 General Modulation Equations for Propagating 
Droplets 

With the results of the previous sections, we now have developed sufficient 
tools to derive the droplet modulation equations. Previous attempts to do 
this have been limited either to a partial set of equations for V and u only, 
computed numerically Hoefer et al. 2012 , or stationary droplet equations 
Bookman and Hoefer| |2013 . The results of Sections [2] and [3] enable us 


to determine the slow evolution of all six soliton parameters due to the 
perturbation p in eq. ([!]). The calculation, not presented, requires some care 
in preserving the appropriate asymptotic relations (14). From expressions 
([5]), (12), and (13), we obtain the droplet soliton modulation equations 


= h />' ' ,)sech (" ^ l) ped *+y 


■Vo; 

^-7 + 2 ^ 


sech ( p -) ppedx, 


sech ( p -) p$dx, 


(23) 

(24) 


6j = — — I sech [ p -) pedx, 


V =- 


U! 

7r 


u 


-V- 

2 


(Vy) 

poj 


Cp ) sech ( p - ^ j p©dx - 7- 


(25) 


sech p - 

ui 

(26) 


where the over dot denotes differentiation with respect to T. This general 
set of equations is the main result of this work. They are asymptotically 
valid when 


0<e<l, T<Ce 


-1 


0 < IVI « 00 « 1. 


(27) 


The perturbation components p@, p$ are to be evaluated with the approx¬ 
imate droplet solution (12), (13). Some explanation of the perturbation 


components and the unit vectors p, Cp is warranted. The magnetization m 
has the unit sphere § 2 as its range. We can therefore define the standard, 
right-handed, orthonormal, spherical basis |r,4»,©| for M 3 where 


r = m 


pp$dx, 
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is the radial unit vector, <l> is the azimuthal unit vector, and 0 is the polar 
unit vector. The components of p in this “magnetization centered” basis are 


0 = ? • p, p# = $ ■ p, p© = © ■ p. 


(28) 


On the other hand, the domain R 2 has the standard orthonormal, polar basis 
{/>, ip} where p and <p are the radial and azimuthal unit vectors, respectively. 
This corresponds to the “domain centered” basis. It is important not to 
confuse the domain R 2 and range § 2 of m. 

In this general formulation, we have neglected spatial inhonrogeneity 
of the perpendicular magnetic field magnitude ho- In Section 5 5.1 


we 


will incorporate inhomogeneity as a perturbation with nonzero p<j> compo¬ 
nent. Even with an inhomogeneous magnetic field, the total frequency of the 
droplet Q(T) is 


Q(T) = ho(e£(T),T) + w(T) + e6 0 (T). 


(29) 


We see that variations in the initial phase <I>o provide a higher order cor¬ 
rection to the droplet frequency. Additionally, the second term on the right 


hand side of eq. (24) is a higher order correction to the droplet’s total ve¬ 
locity These higher order corrections have proven to be of fundamental 
importance in the study of stationary droplets Bookman and Hoefer |2013 
and beyond, see, e.g., Ablowitz et al. [2011 for an application to NLS dark 
solitons. 

While quite general, these equations do not treat all perturbations. It 
is important to note that the solvability condition which gives rise to these 
equations applies for those perturbations whose temporal dependence is on 
a slow time scale. In the sections that follow, we consider a range of physical 
perturbations that meet this criterion in order to demonstrate the versatility 
of this approach. However, some physical scenarios (such as an applied field 
varying rapidly in time) might not satisfy this assumption. Such pertur¬ 
bations may be regular perturbations and not induce dynamics within the 
family of solitons so they will not be further discussed. 


5 Applications To Perturbed Systems 

In this section we analyze a range of perturbations to demonstrate the ver¬ 
satility of this framework as well as to provide physical insights into droplet 
dynamics. 

5.1 Slowly Varying Applied Field 

In practical applications, the magnetic field will typically have some spatial 
variation whose scale is much larger than the scale of the droplet, i.e., the 
exchange length divided by u ;. For this, we assume that /iq = ho(ex, et), 0 < 
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e/u <C 1. This inhomogeneity is best treated by introducing an appropriate 
perturbation p in eq. ([!]). Expanding ho about the soliton center, £, 

h 0 (ex, et) = h 0 (e£, et) + eV/i 0 | x=€ • (x - £) + O (e 2 ) , (30) 


where V represents the gradient with respect to the slow variable X = ex. 
Inserting the expansion (30) into the cross product —m X (hoz) from eq. ([!]) 
introduces the perturbation 


Pe = 0 and p$> = ( Vh o • p. 


(31) 


Substituting these into eqs. {23)-(26) leads to Newton’s second law for the 
droplet center 


<?£ _ (IV 
~dP “ 


—wV/io, 


(32) 


Note that V here represents the gradient with respect to the fast variable x, 
distinguishing it from V. The phase d>o and frequency u are unchanged by 
the field gradient. 

A favorable comparison of direct numerical simulations for eq. ([I]) (see 
Appendix [A]) with the solution to (32) is shown in Fig. [2] We now demon- 



t 


Figure 2: Acceleration of the droplet due to the inhomogeneous magnetic 
field ho = 0.5 — 10 -4 x with w(0) = 0.1 and | V(0)| = 0. The exact solution 
to eq. ( |32| ) (solid) compares favorably to direct numerical simulations of the 
PDE (dashed). 


strate that the explicit equation (32) agrees with the previous result in Hoefer 
et al. |2012j obtained by perturbing conservation laws and integrating the 
equations numerically. Previously, the nontrivial dynamical equation was 
^ = —ATVho. We can transform this equation into eq. (32) by using the 
explicit formulae (18), (16) for A f and V. Since 4^ = 0 and J\f depends only 


on cj, = 0. Then 


dP _ AA 3 / 2 dX _ m e s d 2 £ _ A/" 

dT ~ df ~ ~ e h ° 


(33) 
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This is exactly (32). The particle-like droplet with mass m e g- in eq. (22) 


experiences a conservative force due to the potential J\Tho. This interpre¬ 
tation is consistent with the analysis of the effective mass derived from the 
kinetic energy in Section [3] Furthermore, it demonstrates that a droplet in a 
magnetic held gradient behaves effectively like a single magnetic dipole with 
net dipole moment J\T. 

The effect of an inhomogeneous magnetic held on a massive two-dimensional 
droplet is markedly different from its effect on a one-dimensional droplet 
|l998] and a vortex Papanicolaou and Tomaras [11991 j. A 


Kosevich et al. 


one-dimensional droplet experiences periodic, Bloch-type oscillations for a 
magnetic held with constant gradient, while a magnetic vortex exhibits mo¬ 
tion perpendicular to the held gradient direction. 


5.2 Damping 


In Hoefer et al. 12012 , it was observed that the droplet accelerates as it 


collapses in the presence of damping alone. The framework presented here 
offers an analytical tool to understand this slightly counterintuitive result, 
namely that damping can cause the otherwise steady droplet to speed up. 
The relevant contributions to eq. [l] are 


pe = —(oj + ho — V • V<h) sin(0) and p<j> = —V • V0 


(34) 


where the small parameter e is the Landau-Lifshitz magnetic damping pa¬ 
rameter, usually denoted a. In many practical situations, the damping pa¬ 
rameter is quite small. 


Evaluation of equations (23)-(26) with these perturbations yields two 
nontrivial equations 


du o / i \ 

— = u (w + ho) 

dV 

= wV (to + 2ho) ■ 


(35) 

(36) 


These equations are again consistent with the numerical, perturbed conser¬ 
vation law approach taken in Hoefer et al. [2012 when evaluated at the 
approximate solution. 

We observe that the right hand sides of the modulation equations are 
both positive for ho > — uj/2. Hence, the frequency and velocity increase. 


Equation (35) can be interpreted as a dynamical equation for the droplet’s 


mass m e ff (eq. (22)). The mass is decreasing at a faster rate than the velocity. 
In light of the interpretation given in Section ^5.1 even though the droplet 
is losing energy, it sheds mass fast enough that its acceleration is not a 
contradiction. In Fig. [3] we see good agreement between the modulation 
theory and full micromagnetic simulations. 


15 


































t t 


Figure 3: The evolution of droplet frequency (a) and velocity (b) due to 
damping for both numerical solutions of eqs. ( |35| ), ( |36| ) (solid) and direct 
numerical simulations of eq. ([!]) (dashed) when e = a = 0.01, Hq = 0.5, 
w(0) = 0.1 and |V(0)| = 0.01. 


Since Eq. (35) decouples in this system, an analytical solution can be 


found. Elementary application of partial fractions yields an explicit solution 
in terms of the Lambert W-function; however, the analysis is significantly 


simplified when ho = 0. In this case, the analytical solution to Eqs. (35)-(36) 


is 


u(t) = 

V(t) = 


wo 


y/l — 2acugt 

Vo 

a /1 — 2auQt ’ 


(37) 

(38) 


where Wo is the initial precession frequency and Vo the initial velocity. These 
expressions reveal two facts: a clear time of breakdown for modulation the¬ 
ory and the existence of an adiabatic invariant. Dividing Eq. (37) by the 


components of Eq. (38) demonstrates that the quantities u/V x and U)/V y are 
constant in time. 


5.3 NC-STO and Spatially Inhomogeneous Applied Field 

So far, the examples we have chosen to focus on have not included higher 
order contributions via the phase $0 and the second term of the equation 
for the droplet center £. But many perturbations and physical behaviors 
cannot be investigated without these higher order terms. Consider the more 
complex system of a nanocontact spin-torque oscillator (NC-STO), in which 
a polarized spin current exerts a torque on the magnetization, the spin trans¬ 
fer torque Berger |1996 , Slonczewski |1996 . This forcing can be confined to 
a localized region via a nanocontact Tsoi et al. |1998 , Slonczewski |1999| , 
Rippard et al. |2004 . Perturbations of this sort lead to dynamics within all 
the parameters of the droplet. In addition to spin torque, a droplet in a 
NC-STO also experiences damping and it is precisely the balance between 
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the two that leads to the stable droplet observed in experiments. Here, we 
consider the addition of weak spatial inhomogeneity of the applied magnetic 
field. For simplicity we restrict onr consideration to a constant magnetic 
field gradient. 

This investigation has broader implications for the practical use and un¬ 
derstanding of droplets in real devices. We show in this section that these 
three physical effects influence the system in competing ways, which can 
balance, allowing for the existence of stable droplets. Alternatively, a strong 
enough field gradient can push the droplet out of the NC-STO, giving rise 
to a previously unexplained drift instability Hoefer et al. 120101. As seen in 
Sec. ^3.2 damping decreases the effective mass of the droplet. In Sec. ^3.1 


it was shown that a field inhomogeneity accelerates the droplet while leaving 
the mass of the droplet unaffected. The inclusion of forcing due to spin trans¬ 
fer torque in a nanocontact opposes both of these effects. The spin torque 
increases the droplet mass and generates an effective restoring force that 
centers the droplet in the nanocontact region Bookman and Hoefer [2013], 
Hence, there can exist a delicate balance between all of these effects: the NC- 
STO restoring force balancing the potential force due to the field gradient 
and the mass loss due to damping balancing the mass gain due to spin- 
torque. Previous studies have been unable to identify when such a balance 
occurs and when it fails. Here, we analytically demonstrate stable droplets 
as fixed points of the modulation equations with all of these perturbations. 

Because the perturbation components p© and p<j, appear linearly in the 
modulation equations (p3|)-([26|), we can simply add the field inhomogeneity 


eq. (31) and damping eq. (|34|) perturbations to those due to spin torque Stiles 


and Miltat [2006]. Due to the presence of three different perturbations, we no 
longer scale the perturbation p in eq. ([Tj) by the single parameter e. Rather, 
we set e = 1 and introduce the small parameters in p© and p<j> directly. The 
perturbation components are 

p© = —a {uj + ho — V • V4>) sin 0 + a — r ) sin 0, (39) 

p<j> = (V/io ■ p) p — aV • V0. (40) 

The nanocontact where spin torque is active is assumed to be a circle with 
radius p*. The coordinate r in the argument of the Heaviside function Jf 1 
is measured from the center of the nanocontact, which differs from the co¬ 
ordinates p and p which are measured from the center of the droplet. For 
simplicity, we have neglected the spin torque asymmetry that introduces an¬ 
other parameter into the analysis but does not appear to have a significant 


on the dynamics 

Hoefer et al. 

|2010 . Experiments 

Ylohseni et al. 

Macia et al. 

2011 and analysis 

Hoefer et al. 

(2010 , 

Bookman and 


cr (proportional to current), is roughly order 1 for the existence of droplets 
to be satisfied. Thus 0 < a ~ a <C 1. The magnetic field is assumed to be 
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linear 


ho = a + bx , |6|<Cu;. (41) 

We can restrict to droplet motion in the x direction only. Insertion of the 


perturbations in Eqs. (39), (40) into the modulation equations (23)- (26) 


results in the following system 
abV aV 


4>n = 


—-— / cos(</?)sech 2 p -dx 


2oo 2 47t 




1 


LO 


i = V — — + — [ cos(y?)sech 2 ( p — — ) dx 


u) 


2 vr 


M<P* 


00 


Co = aio 2 (lo + a) — 


GOO 
47T 


V = —boo + aV lo (lo + 2a) — 


|x|<p* 

aVoo 


sech z ( p -) dx 


LO 


(3 poo + cos(2(^) — 1) 


47T 


'|x|<p* 


(42) 

(43) 

(44) 


sech p - 

00 


(45) 


where we set £ = V = V x , and the over dot denotes differentiation 
with respect to t. None of the right hand sides in the equations above 
depend explicitly on the parameter 4>o so that the dynamics of the remaining 
parameters can be considered separately. We ignore the evolution of 4>o for 
the remainder of the analysis noting that 4>o corresponds to a small frequency 


shift as in Eq. (29) that can be obtained from the evolution of the other 


parameters by insertion into Eq. (42). 


There is a complex interplay between the many small parameters in this 
problem. Since we do not have access to an exact analytical solution, it 
is necessary that these perturbations dominate over the error terms in our 
approximate solution, while still remaining small. Since we have \V\ < lo 2 
to keep an overall consistent error estimate for the approximate droplet, we 
require that a, a <C lo. The variation in the applied held, b, is a more subtle 
and the appropriate scaling will be determined by directly computing fixed 
points. 

The stationary droplet without a held gradient is stable when centered 

20101 , 


on the nanocontact Hoefer et al. 


Bookman and Hoefer 2013 . This 


results from an analysis of the stationary modulation equations which exhibit 
an attractive, stable hxed point. Taking ^ = 6 = 0, the modulation 
equations (|43[)-(|45|) for propagating droplets exhibit the same hxed point 
(£,lo,V) = (0,cu*,0) when damping balances forcing corresponding to the 
current 


a 


a 


2 (a T lo *) 


1 + 00 ] 


^log ^ sech 


P* - W ) ) + P* tanh ( P* ~ Jr 


UJ* 


(46) 


dx 
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There is a saddle node bifurcation as a is increased, with cn* = w*(cr) cor¬ 
responding to the stable branch. For o sufficiently large, the stable branch 
quickly approaches 


w* = p* 1 +arctanh P * 2 +0 (p* 3 ) , p* » 1, 0 < w*—p* 1 <C 1. 

(47) 

Near the critical value a = 2 aa, where the second term is small, the asymp¬ 
totic form is 


w* = 


P * 1+ ~ P* 2 + (jjr + ln 2 ) P* 3 +° {p* 4 )> ~^~ 1 =° ( P * x ) • 


(48) 

Linearizing equations (43)- (45) about this fixed point, the Jacobian matrix 
is given by 


J(0,cj*, 0) 

Ai 

A 2 

A3 


/ Ai 0 1 \ 

0 a 2 0 , 

V 0 0 As / 


1 


—-crp*u;*sech p*-, 


1 




. 1 

—craw* + Ai + -ffw, 
—2 aw* + A 2 — Ai- 


^tanh 



(49) 

(50) 

(51) 

(52) 



fixed point is stable. The critical forcing value a = aa, below which the 
droplet may be unstable could be considered as an estimate for the minimum 
sustaining current of a droplet Hoefer et al. [2010|. Note, however, that this 


is a dubious estimate due to a;* — p* 1 not being a small quantity. Utilizing 


the approximation from Eq. (48), we find 


Ai - -- + O (p* 2 ) , A 2 ~ Ai, A 3 - + ^ _cr ) P * 2 + ® 4 ) ’ 


2 aa 


- 1 


a 


(53) 


We now turn our attention to the case of a small field gradient 0 < |6| -C 
1, where we observe the persistence of the droplet fixed point for very small 
16|. These fixed points exist as a balance between the expulsive force provided 
by the field gradient and the attractive force provided by the nanocontact. 
This attraction manifests in the evolution of £ and so this balance can also 
be viewed as a balance between leading order effects (in V) and higher order 
effects (in £). Unlike the stationary fixed point, exact analytical expressions 


= O (ft 1 ). 
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Figure 4: Fixed points from modulation theory, exact (solid) and approxi¬ 
mate Eq. (55) (dashed), and direct numerical simulation of Eq. ([T]) (circles) 
when a = a = 0.01, a = 0.5, /?* = 12. In this case, the parameter V 
cannot be extracted from direct numerical simulations without additional 
assumptions (see Appendix 0 Accordingly, this data is not presented in 
(c). 


for the fixed point cannot be found since the droplet is no longer centered on 
the nanocontact (£ ^ 0). Nevertheless, we can obtain the approximate form 


for these fixed points as follows. The structure of J in Eq. (49) yields very 


simple predictions in the regime of small field gradient. The key observation 


here is that the system of Eqs. (43)-(45) can be written as 


w I =F{£,u,V)-b I 0 

kV ) \Uj 


(54) 


By virtue of the stationary fixed point, F satisfies F(0, w*, 0) = 0. We now 
seek a fixed point that slightly deviates from the stationary one according 
to £ = + • • •, u = w* + buii + • • • and V = bV\ + ■ ■ ■. Expanding and 


equating the right hand side of Eq. (54) to zero gives the correction 



= J(0,w*,0) 


-i 



Aio;* A1A3 


CJ* 

A 3 



(55) 


— 2 Q+cr In 2 


where the approximations (48) and (53) were used to obtain the large p * 
estimate. These approximations are valid so long as a is more than p~ 2 
away from the critical value 2a/ In 2 ~ 2.9a. Otherwise, higher order terms 


in Eq. (53) would need to be considered. 


As summarized in Fig. [l] these simple expressions make predictions in 
good agreement with the fixed points found by numerical continuation in b 
and those observed in long time micromagnetic simulations of Eq. Q with 
perturbations (39) and (40). The Jacobian matrix of Eqs. (43)-M5|) can also 


be numerically evaluated, showing that all eigenvalues are negative, until 
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continuation breaks down when one eigenvalue reaches zero. After this bi¬ 
furcation, we do not find any fixed points. The condition of this eigenvalue 
reaching zero then corresponds exactly to the crossover where the attractive 
nanocontact is no longer strong enough to balance the expulsive force sup¬ 
plied by the field gradient. If we assume that a field gradient is strong enough 
to move the droplet an order one distance, still small relative to the droplet 
radius cj" 1 , then we obtain a typical field gradient scaling b ~ aa /2p*. This 
field gradient is very small. For the example studied here, b ~ 10 -6 com¬ 
pared to the NC-STO forcing magnitude a = 10 -2 . This demonstrates that 
droplet attraction due to spin torque is weak relative to droplet acceleration 
due to field inhomogeneity. A strong enough field gradient, on the scale 
of aa/2p *, can eject the droplet from the nanocontact, causing a drift in¬ 
stability previously observed in numerical simulations Hoefer et al. [20101. 
Additionally, the associated velocity scale from Eq. (55) is V = bV\ ~ a /2 


which is much smaller than oj as required for this order of accuracy of the 
approximate droplet. 


6 Interacting Droplets 

An intriguing, indeed defining, aspect of solitary wave dynamics is their in¬ 
teraction behavior. Reference Maiden et al. 2014 undertook a numerical 


investigation of two interacting droplets by varying droplet parameters and 
quantifying the properties of the solution post-interaction. It was found 
that the relative phase difference between the two droplets plays a funda¬ 
mental role, controlling whether the interaction is attractive or repulsive. 
The attractive interactions studied were strongly nonlinear, hence a pertur¬ 
bation theory would be insufficient to study the full complement of observed 
phenomena. Nevertheless, we can gain insight into the nature of the interac¬ 
tion (attractive/repulsive) by studying two well-separated droplets perturba- 
tively, with the small parameter being the inverse of the droplet separation. 
This approach is well-known and has been applied successfully to, for exam¬ 


ple, NLS-type models Zhu and Yang 2007 , Ablowitz et al. 2009 


In full generality, the perturbations arising from this analysis are com¬ 
plex. However, since the validity of these equations is strongly dependent on 
the separation of the two droplets, we only expect these equations to be valid 
over short time scales. Hence, we only seek to describe the initial behavior of 
two stationary, weakly overlapping droplets. As the interaction immediately 
excites propagation of the two droplets, this will not model the behavior for 
t > 0. Nevertheless, these assumptions make it possible to describe much 
of the behavior observed in full numerical simulations Maiden et al. |2014|. 
The initial configuration places one droplet on the left (subscripted 1) and 
another droplet (subscripted 2) a distance d away along the x— axis. We 
define the relative phase difference A4> = 4>2 — 3>i, which will emerge as 
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an important quantity in the modulation equations. Considering the mod¬ 
ulation equations for two weakly interacting droplets with motion in the x 
direction at the initial time only yields 


UJ 


$o,fc = ~ 7 T cos(AT) / /C fc (x)dx, 

27T . W2 


(56) 


ik = —(—l) fe+1 sin(AT) f /Cfc(x)sech ( p -) cos <pdx, (57) 

2n J r 2 V uJ 

hi, = — — (—l) fc+1 sin(AT) f /Cfc(x)sech (p -^ dx, (58) 

47T J R 2 V w / 

/.,3 r 

(59) 


14 =—cos(AT) / /Cfc(x) cos <pdx, 

7T ./l»2 


where 


/Cfc(x) = sech -^ sech ^p 


w 


2secfr ( p-J —— tanh ^p- 

(60) 


The integration kernel ICk depends on the separation between the two droplets 
through pk = \/(x + (—l) fc d) 2 + y 2 . Utilizing this framework, we can now 
offer some insight into the nature of two interacting droplets. 

6.1 Attraction and Repulsion 

The attractive or repulsive nature of two droplets can be understood by 


considering Eq. 59 As AT varies, the sign of cos(AT) is clear. Determining 


the initial direction of motion, right or left, of the droplet comes down to 
determining the sign of the integral term in (59). Figure [5} left shows the 



0.05 0.1 0.15 0.2 0.25 

u>0 



0.05 0.1 0.15 0.2 0.25 

OJQ 


Figure 5: Left: Initial acceleration for varied initial uj and several values of 
separation. Right: Plot of 6j\ as a function of initial ujq for several values of 
separation. In both, the initial relative phase was AT = 1. 

numerical evaluation of the right hand side of V\ (droplet on left) when AT = 
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1 < vr/2, leading to positive values only. Thus, the left droplet experiences 
a positive acceleration to the right, towards the other droplet when |A3?| < 
7 t/ 2. Since the kernel exhibits the symmetry with respect to droplet choice 


/Ci (x,y) = /c 2 (— x,y), the integral in (59) for the right droplet, k = 2, has 


the opposite sign. The right droplet experiences a negative acceleration to 
the left when |A3>| < ir/2. Therefore, two droplets are attractive when 
| A3?| < vr/2, i.e., when they are sufficiently in phase. Similarly, when vr/2 < 

| A3? | < 7r, the signs of V/ are reversed and the droplets move away from 
each other. Thus, two droplets are repulsive when they are sufficiently out 
of phase. 

As was noted in Maiden et al. 12014] by a nonlinear method of images, the 
attractive or repulsive nature of two droplets with the special initial values 
A3> = 0 or AT = 7r describes the behavior of a single droplet near a mag¬ 
netic boundary with either a free spin (Neumann type) boundary condition 
or a fixed spin (Dirichlet type) boundary condition, respectively. The anal¬ 
ysis presented here confirms this fact for any droplet that weakly interacts 
with a magnetic boundary. Such behavior was observed in micromagnetic 


simulations of a droplet in a NC-STO, nanowire geometry Iacocca et al. 


12014 


6.2 Asymmetry 

Despite a highly symmetric initial condition, an asymmetry was observed in 
so-called “head-on collisions" of two droplets in Maiden et al. [20141. The 
frequency equation ( |58[ ) provides an explanation of this in the limit of very 
small velocities. Figure |5j right contains the relevant information. Since 
w\ appears to be negative always when sin(AT) > 0. In the numerical 


experiments in Maiden et al. 2014 were done over the range A3? = 0 to 


A3> = 7r, this was always the case. Again using that ICi(x,y) = IC 2 (—x,y), 
it can be seen that the integrals involved in computing cLq and 002 are equal. 
Hence the sign of Uk is determined by (—l) fc+1 , and the signs of Cj 1 and 
UJ 2 w hl always be opposite. For the parameters discussed here, this means 
that the frequency decreases for the droplet on the left and increases on the 
right. This change in droplet structure is asymmetric because a reduced 
(increased) frequency implies larger (smaller) droplet mass and corresponds 
precisely with the observations of Maiden et al. 12014 . 


6.3 Acceleration 


The discussion of attraction and repulsion in Section [t|p.l| suggests that the 
boundary between the two behaviors is A3> = n/2. But this does not agree 
with numerical experiments where the crossover A3> was found to vary with 
the initial droplet parameters Maiden et al. [201 41. To offer an explanation 
for this, we consider the total acceleration of the initial droplets, i.e., £*.. 
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This incorporates higher order information not included in Vj,. Since the 
full modulation equations for interacting droplets when V 7 ^ 0 are complex, 
we do not examine for all values of AT. However, at AT = we know 
Vk = To = 0, (since cos(AT) = 0) and those terms will not contribute, which 
simplifies the calculation. Figure [ 6 ] shows the initial, total droplet acceler¬ 
ation £ 1 , evaluated numerically, as the initial frequency and separation are 
varied. The variable sign of this quantity as parameters change demonstrates 
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Figure 6 : Numerical evaluation of £1 initially for AT = vr/2, variable droplet 
separation d and frequency ujq. There is not one sign of acceleration, i.e., 
the left droplet can be repelled or attracted to the right droplet depending 
on the choice of parameters. 

that subtle, higher order effects causes the crossover value of AT to deviate 
from its nominal value n/2. 

7 Conclusion 

The primary contribution of this work is a general framework for investigat¬ 
ing perturbations of droplet solitons. Actual physical devices used to create 
and manipulate droplets are quite complex, incorporating a number of phys¬ 
ical effects. Therefore, having a tractable, analytical theory to describe both 
the motion and precession of droplets due to physical perturbations is quite 
valuable. The examples presented here are meant to demonstrate the ver¬ 
satility and power of this tool. Additionally, the application of this theory 
to the NC-STO provides several insights into the behavior of experimentally 
observed dissipative droplets. In particular, the dissipative droplet is shown 
to be robust in the presence of weak held gradients, but can be ejected from 
the nanocontact if the held gradient is too large, providing an explanation 
for a previously observed drift instability. These observations open possi¬ 
ble mechanisms for generating a current of solitons which could serve as a 
mechanism for information transfer. 

As demonstrated by examples in the preceding sections, many perturba¬ 
tions excite evolution of higher order parameters (overall phase and position). 
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This subtle information proves to be of fundamental importance for several 
perturbations considered. As shown by our derivation of the modulation 
equations for a general class of Hamiltonian systems, the higher order pa¬ 
rameter dynamics emerge when the generalized nullspace of the linearized 
evolution operator is incorporated. Due to the existence of an approximate 
analytical form for the propagating droplet, we are able to completely char¬ 
acterize this nullspace and hence recover the droplet modulation equations 
in a convenient form. 

A number of physical perturbations can now be investigated within this 
framework. Future developments of the modulation theory itself could be 
performed in the context of a different family of magnetic droplet solutions. 
One example is the weakly nonlinear droplet Ivanov et al. [2001], Recent mi- 
cronragnetic simulations have found rotating and processing localized waves 
in NC-STOs with non-trivial magnetostatic contributions Finocchio et al. 


[)2013 . The invariance of equation ([I]) when e = 0 with respect to rota¬ 
tion of the domain and an analysis of conserved quantities Papanicolaou and 
Tomaras [1991], suggests that there may be rotating and processing solitary 
wave solutions. The modulation theory developed in this work could be 
extended to such solitary wave solutions. 


Acknowledgment 

The authors gratefully acknowledge support through an NSF CAREER 
grant. 


References 

M. J. Ablowitz and H. Segur. Solitons and the inverse scattering transform, 
volume 4. SIAM, 1981. 

M. J. Ablowitz, T. P. Horikis, S. D. Nixon, and Y. Zhu. Asymptotic analysis 
of pulse dynamics in mode-locked lasers. Stud. Appl. Math., 122(4) :411- 
425, 2009. 

M. J. Ablowitz, S. D. Nixon, T. P. Horikis, and D. J. Frantzeskakis. Per¬ 
turbations of dark solitons. Proc. R. Soc. A-Math. Phys., 467(2133):2597 
-2621, 2011. 

I. M. Babich and A. M. Kosevich. Relaxation of bloch oscillations of a 
magnetic soliton in a nonunifornr magnetic held. Low Temp. Phys., 27: 
35-39, 2001. 

L. Berger. Emission of spin waves by a magnetic multilayer traversed by a 
current. Phys. Rev. B, 54(13):9353, 1996. 


25 




















L. D. Bookman and M. A. Hoefer. Analytical theory of modulated magnetic 
solitons. Phys. Rev. B, 88(18): 184401, 2013. 

S. Chung, S. M. Mohseni, S. R. Sani, E. Iacocca, R. K. Dumas, T. N. Anh 
Nguyen, Ye Pogoryelov, P. Iv. Muduli, A. Eklund, M. Hoefer, and J. Ak- 
erman. Spin transfer torque generated magnetic droplet solitons (invited). 
J. Appl. Phys., 115:172612, 2014. 

A. Fert, V. Cros, and J. Sampaio. Skyrmions on the track. Nature nanotech¬ 
nology, 8(3): 152—156, 2013. 

G. Finocchio, V. Puliafito, S. Komineas, L. Torres, O. Ozatay, T. Hauet, and 
B. Azzerboni. Nanoscale spintronic oscillators based on the excitation of 
confined soliton modes. J. Appl. Phys., 114(16): 163908, 2013. 

M. Grillakis, J. Shatah, and W. Strauss. Stability theory of solitary waves 
in the presence of symmetry, ii. J. Fund. Anal., 94(2):308-348, 1990. 

M. A. Hoefer and M. Sommacal. Propagating two-dimensional magnetic 
droplets. Physica D, 241(9):890-901, 2012. 

M. A. Hoefer, T. J. Silva, and M. W. Keller. Theory for a dissipative droplet 
soliton excited by a spin torque nanocontact. Phys. Rev. B, 82(5):054432, 
2010 . 

M. A. Hoefer, M. Sommacal, and T. J. Silva. Propagation and control of 
nanoscale magnetic-droplet solitons. Phys. Rev. B, 85(21):214433, 2012. 

E. Iacocca, R. K. Dumas, L. Bookman, M. Mohseni, S. Chung, M. A. Hoe¬ 
fer, and J. Akerman. Confined dissipative droplet solitons in spin-valve 
nanowires with perpendicular magnetic anisotropy. Phys. Rev. Lett., 112 
(4):047201, 2014. 

B. A. Ivanov and V. A. Stephanovich. Two-dimensional soliton dynamics in 
ferromagnets. Phys. Lett. A, 141(l):89-94, 1989. 

B. A. Ivanov, C. E. Zaspel, and I. A. Yastremsky. Small-amplitude mobile 
solitons in the two-dimensional ferromagnet. Phys. Rev. B, 63(13):134413, 
2001. 

J. P. Keener and D. W. McLaughlin. Solitons under perturbations. Phys. 
Rev. A, 16:777, 1977. 

P. G. Kevrekidis, D. J. Frantzeskakis, and R. Carretero-Gonzalez. Emergent 
nonlinear phenomena in Bose-Einstein condensates : theory and experi¬ 
ment. Springer, 2008. 

Y. S. Kivshar and B. A. Malomed. Dynamics of solitons in nearly integrable 
systems. Rev. Mod. Phys., 61(4):763, 1989. 


26 



A. M. Kosevich, B. A. Ivanov, and A. S. Kovalev. Sov. Sci. Rev., Sect. A, 
Phys. Rev., 6:171, 1986. 

A. M. Kosevich, B. A. Ivanov, and A. S. Kovalev. Magnetic solitons. Physics 
Reports, 194(3): 117-238, 1990. 

A. M. Kosevich, V. V. Gann, A. I. Zhukov, and V. P. Voronov. Magnetic 
soliton motion in a nonuniform magnetic field. J. Exp. Theor. Phys., 87: 
401-407, 1998. 

F. Macia, D. Backes, and A. D. Kent. Stable magnetic droplet solitons in 
spin-transfer nanocontacts. Nat. Nanotechnol., 9(12):992-996, 2014. 

M. D. Maiden, L. D. Bookman, and M. A. Hoefer. Attraction, merger, 
reflection, and annihilation in magnetic droplet soliton scattering. Phys. 
Rev. B, 89(18):180409, 2014. 

S. M. Mohseni, S. R. Sani, J. Persson, T. N. Anh Nguyen, S. Chung, Y. Pogo- 
ryelov, P. K. Muduli, E. Iacocca, A. Eklund, R. K. Dumas, S. Bonetti, 
Deac A., M. A. Hoefer, and J. Akerman. Spin torque-generated magnetic 
droplet solitons. Science, 339(6125) :1295-1298, 2013. 

S. M. Mohseni, S. R. Sani, R. K. Dumas, J. Persson, T. N. Anh Nguyen, 
S. Chung, Ye. Pogoryelov, P. K. Muduli, E. Iacocca, A. Eklund, and J. Ak¬ 
erman. Magnetic droplet solitons in orthogonal nano-contact spin torque 
oscillators. Physica B, 435:84-87, 2014. 

N. Nagaosa and Y. Tokura. Topological properties and dynamics of magnetic 
skyrmions. Nat. Nanotechnol., 8(12):899-911, 2013. 

N. Papanicolaou and T. N. Tomaras. Dynamics of magnetic vortices. Nucl. 
Phys. B, 360(2):425-462, 1991. 

B. Piette and W. J. Zakrzewski. Localized solutions in a two-dimensional 
landau-lifshitz model. Physica D, 119(3):314-326, 1998. 

J. F. Pulecio, P. Warnicke, S. D. Pollard, D. A. Arena, and Y. Zhu. Co¬ 
herence and modality of driven interlayer-coupled magnetic vortices. Nat. 
Commun., 5, 2014. 

W. H Rippard, M. R Pufall, S. Kaka, S. E Russek, and T. J Silva. Direct- 
current induced dynamics in Co 9 oFeio/NigoFe 2 o- Phys. Rev. Lett., 92(2): 
027201, 2004. 

A. Sanchez and Bishop A. R. Collective coordinates and length-scale compe¬ 
tition in spatially inhomogeneous soliton-bearing equations. SIAM Rev., 
40(3):579-615, 1998. 


27 



J. C. Slonczewski. Current-driven excitation of magnetic multilayers. J. 
Magn. Magn. Mater., 159(1):L1~L7, 1996. 

J. C. Slonczewski. Excitation of spin waves by an electric current. J. Magn. 
Magn. Mater., 195:L261-L268, 1999. 

M. D. Stiles and J. Miltat. Spin transfer torque and dynamics. In B. Hille- 
brands and A. Thiaville, editors, Spin Dynamics in Confined Magnetic 
Structures III. Springer, Berlin, Heidelberg, 2006. 

M. Tsoi, A. G. M. Jansen, J. Bass, W.-C. Chiang, M. Seek, V. Tsoi, and 
P. Wyder. Excitation of a magnetic multilayer by an electric current. 
Phys. Rev. Lett., 80(19):4281-4284, 1998. 

V. Uhlif, M. Urbanek, L. Hladik, J. Spousta, M.-Y. Inr, P. Fischer, N. Eibagi, 
J. J. Kan, E. E. Fullerton, and T. Sikola. Dynamic switching of the spin 
circulation in tapered magnetic nanodisks. Nat. Nanotechnol., 8(5):341- 
346, 2013. 

N. G. Vakhitov and A. A. Kolokolov. Stationary solutions of the wave equa¬ 
tion in a medium with nonlinearity saturation. Radiophys. Quantum Elec¬ 
tron., 16(7):783-789, 1973. 

M. I. Weinstein. Modulational stability of ground states of nonlinear 
schrodinger equations. SIAM J. Math. Anal, 16(3):472-491, 1985. 

J. Yang. Nonlinear Waves in Integrable and Nonintegrable Systems. Society 
for Industrial and Applied Mathematics, 2010. 

X. Z. Yu, Y. Tokunaga, Y. Kaneko, W. Z. Zhang, K. Kimoto, Y. Matsui, 
Y. Taguchi, and Y. Tokura. Biskyrmion states and their current-driven 
motion in a layered manganite. Nat. Commun., 5, 2014. 

Y. Zhu and J. Yang. Universal fractal structures in the weak interaction of 
solitary waves in generalized nonlinear schrodinger equations. Phys. Rev. 
E, 75(3):036605, 2007. 


A Numerical Method 

The numerical simulations (micronragnetics) we conducted incorporated a 
periodic Fourier psuedospectral spatial discretization. Unless otherwise stated, 
the spatial domain was [—50, 50] X [—50, 50], sufficiently large so that the per¬ 
turbed solitary waves were well-localized within it. In each spatial dimen¬ 
sion, 2 9 grid points were used. The time-stepping was done using a version 
of the Runge-Kutta-Fehlberg algorithm, modified so that the magnetization 
maintained unit length at every grid point and each time step. 
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The velocity, V, was extracted from numerical data by computing the 
center of mass, £(t) = f K2 x(l — m$(x, t))dx./Af. We then estimated V = £, 
approximated using a forward difference of £ (t ). This method does not work 
for perturbations which excite higher order changes in £ and we do not 
estimate V in such cases. 

For the precessional frequency oj, the phase of the in-plane magnetization 
(m x ,m y ) was extracted at a point a fixed distance from the center of mass 


Differentiating this phase with respect to time yields Cl(t) as in Eq. (29). 


The precessional frequency, oj, was obtained by subtracting /iq and $o- The 


contribution from <f>o was estimated via the modulation equation (23). An 


alternative method based on computing the conserved quantities in Eqs. (15) 


(16) was used for comparison. The relations for total spin and momentum 


of the approximate droplet Eqs. (18)-(19) were inverted to obtain oj and V. 


The two methods were in good agreement. 


B Approximate Droplet 

As noted in the main body of the manuscript, the derivation of the approx¬ 
imate droplet is significantly simplified by exploiting the invariance of Eq. 
( |11[ ) under rotation of the domain and working in the frame where V y = 0 
and V = V x x = Vx. The derivation proceeds by substituting the ansatz 

0 = © 0 (p) + VQi(p, ip) + O (V 2 ) and $ = $ 0 + VQ^p, p) + O ( V 2 ) (61) 


into Eq. (11). At order O (1), this yields one nontrivial equation, 


d ©o . 1 d ®o , n . (ri \ n 

-w H- 7 -h [OJ - cos) sm(© 0 ) = 0. (62) 

ap z p dp 

This equation has been extensively studied and is known to admit an approxi¬ 
mate solution in the limit of 0 < oj <C 1, ©o = ar ccos (tanh (p — ^))+C ) ( a;2 ) 
Ivanov and Stephanovich 119891, Kosevich et al. |1986 , Bookman and Hoefer 
pl3|. From here on out, we additionally assume that 0<w<l. At order 
O (V), we find 


sin(©o)AT] + ( cos(<p) + 2008 ( 00 )-^-^) —= 0, 


dp J dp 

A©i + (oj cos (0q) — cos (2©o)) ©i = 0. 


(63) 

(64) 


Equation (64) is solved by ©i = 0. Substituting the approximate solution 


for ©o into Eq. (63) 


^A4>i — cos (ip) — 2 tanh ^p — — ^ sec h ^P — — ^ = 0 


(65) 


The residual in Eq. (65) is determined by two considerations. If 1 <C p — oj 


-l 


sech (p — 1 /oj) dominates and the residual is exponentially small. In the 
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other case, i.e. p ~ oo 1 ■ the residual will only be small if A<f>i — cos(^) — 
2tanh (p — 1 /oo) is small since sech (p — l/oo) is O (1). This suggests that 
the boundary condition, linipVd)! = ^x, may be neglected. Assuming d>i 


is separable of the form &i(p,<p) = f(p)cos((p), Eq. (65) simplifies to the 
ordinary differential equation 


!|d + (i-2tanh(p-i))!-i/ = l 


P 




dp 


( 66 ) 


Numerical solutions of Eq. (66) demonstrate that / becomes quite large, 


approximately O (-^) near p = Factoring this into the analysis, we 
change to the coordinate system R = p — and expand / in the series 


f(p) = ^ + — +f2( P ) + 


00 


00 


(67) 


Let L = 4ty7. — 2tanh(ii)d^. Substituting the ansatz in Eq. (67) into Eq. 


dR 2 


(66) yields, 


O^j:L/„ = 0 

° © : = -i 

o(i):Lh= -jk +i+f - +Rd A 


( 68 ) 

(69) 

(70) 


Eq. (68) admits any constant solution. Take /q = A. Substituting this 


expression for /o into Eq. (69), yields L f\ = 0. Thus, any constant solution 


is admissible for f\ as well. Take f\ = B. Substituting these expressions for 


/o and f\ into Eq. (70), yields L /2 = 1 + A. However, constant solutions are 


in the kernel of L, so solvability of this equation requires A = — 1. Similarly, 
solvability at O (oo) requires B = 0. Hence, we take f(p) = —which gives 
rise to the form of the approximate droplet in Eq. ( |13[ ). 

Since enters into Eq. [Tl] only in the form of V<L, the asymptotic ap¬ 
proximations here are valid provided that Vd* remains small. This condition 

re approximation for /, 


is not equivalent <C $o- Given the form of t 

the small gradient condition requires that |V3>| = 
approximation is localized at p ~ w -1 , we require V <C oo. 


V sin<p 


< 1. Since the 
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